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Abstract 

Vicious Brownian motion is a diffusion scaling limit of Fisher's vicious walk model, which is a 
system of Brownian particles in one dimension such that if two of them meet they kill each other. We 
consider the vicious Brownian motion conditioned never to collide with each other, and call it the 
noncolliding Brownian motion. This conditional diffusion process is equivalent to the eigenvalue 
process of a Hermitian-matrix-valued Brownian motion studied by Dyson. Recently O'Connell 
introduced a generalization of the noncolliding Brownian motion by using the eigenfunctions (the 
Whittaker functions) of the quantum Toda lattice in order to analyze a directed polymer model in 
1+1 dimensions. We consider a system of one-dimensional Brownian motions with a long-ranged 
killing term as a generalization of the vicious Brownian motion and construct the O'Connell process 
as a conditional process of the killing Brownian motions to survive forever. 
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I. INTRODUCTION 

Vicious walk model introduced by Fisher [l] is the system of one- dimensional random 
walkers such that, if neighboring walkers meet, they kill each other. Though the model 
looks sinister, what we are interested in is to evaluate the probability that for a finite time- 
interval any neighboring pair of vicious walkers do not meet and thus all walkers survive; in 



other words, the probability that the peace is kept 



-|4|. If we take appropriate continuum 



limit (the diffusion scaling limit), we obtain "vicious Brownian motion" (vicious BM) 0,0]. 
Assume that the number of particles of vicious BM is iV > 2 and write the positions as 
1 < j < N. Then the configuration space of them conditioned never to collide is 

Wjy = {x = (x 1 , x 2l • • • , x N ) G R N : x 1 < x 2 < ■ ■ ■ < x N }, (1) 

which is called the Weyl chamber of type Aat_i in the representation theory [7]. The bound- 
aries of this region dWjy in the iV-dimensional real space IR consists of the hyperplanes 
Xj = Xj + i, 1 < j < N — 1, each of which corresponds to occurrence of collision of the j-th 
and (j + l)-th particles in the vicious BM. If we regard x as a position vector of the N- 
dimensional BM within and dW^ as an absorbing boundary such that when a particle 
hit the boundary it is immediately absorbed, the system is identified with the absorbing BM 
in Wat. The harmonic function, A/iat(x) = ^, =1 cPhx^x) /dx'j = 0, x G Wat, satisfying the 
Dirichlet boundary condition /iat(x) = 0,x G dWjy, is uniquely determined up to a constant 
factor as 

M*0= II ( X k-Xj), xGW^, (2) 

l<j<k<N 

which is equal to the Vandermonde determinant deti^j^^iv^j -1 ]- Then we can show that 
the survival probability up to time t > of iV-particle system of vicious BM starting from 
an initial configuration x6 Wat has the asymptotics 



'iV 



where |x| = x\ + • • ■ + x 2 N (see, for instance, [8|). 

As a matter of course, the survival probability decreases in time for any x G Wjy. The 
result ([3]) implies, however, the decay is slow in the sense that it is not exponential but follows 
the power-law in time. This observation has led us to study the system of BMs conditioned 
never to collide with each other, which we call the noncolliding BM. The fundamental 



properties of this process are the following jg]. (i) Let the transition probability density 
of a single BM be p(t, y\x) = e~( y ~ x ^ I" 11 / V 27ri. Then the transition probability density of 
the absorbing BM in Wat from a configuration x £ Wjv to y G Wjv in time interval t > is 
given by the Karlin-McGregor determinant 9| 



q N (t,y\x.) = det Jp(t,yj\x k )], 

1<J, k<N 



(4) 



or equivalently lOj given by the Harish-Chandra-Izykson-Zuber integral III, ll2|] 

qN (t, y|x) = M*)My) / e-^-^^^du, 

C N JV(N) 



(5) 



where dU is the Haar measure of the space of unitary matrices U(iV) normalized 
1, A x = diag(xi,x 2 , . . . ,x N ), A y = di&g(y 1 ,y 2 , . . . ,y N ), and c N = 



as iu(JV) ^ 



(2 



.N/2' 



Y\.j=i^{j) with the gamma function T(j). (ii) The transition probability density 
of the noncolliding BM, p^(t,y\x), is then given by the harmonic transform of qN{t,y\x.) 



with fl2l) in the sense of Doob 



pjv(t,y|x) 



13| 



My) 

h N (x) 



q N (t,y\x), xjgWjv, £>0. 



(6) 



(iii) If we regard (j6| as a function of t and initial configuration x, it is a solution of the 
following backward Kolmogorov equation 



d_ 

at 



u(t,x) = -Au(t, x) + V log fojv(x) • Vw(t, x) 



1 N d 2 



1 



<9 



l^j^N l<k<N:k^j 



u(t, x) 



(?) 



under the initial condition u(0, x) = <5(x — y) = n,=i <K x j ~ Vj)- 

The direct consequence of the above facts is the following (see, for example, [jjj]). Let 
X(£) = (Xi(t), . . . ,Xjtf(t)) be the iV-particle system of the noncolliding BM. Then it solves 
the system of stochastic differential equations (SDEs) 

dt 



dXj(t) = dBj(t) + J2 



x 3 (t)-x k (ty 



l<j<N, t > 0, 



(8) 



l<k<N:kj±j 

where Bj(t),l < j < N, t > 0, are independent one-dimensional standard BMs. Eq.(|8) 
is nothing but the system of SDEs for Dy son's BM model with the parameter (3 = 2 15 1 



studied in the random matrix theory 



16 



17| . which we will simply call the Dyson model 



in this paper. That is, the nonco 



Hermitian-matrix-valued BM 



10 



liding BM is equivalent to the eigenvalue process of the 

Recently O'Connell introduced a family of diffusion processes of N particles in one di- 
mension, Z M (£) = (Zi(t), . . . , Zjjf(t)), t > with a parameter given by an iV-dimensional 
real vector \x = /i 2 , • • • , Hn) £ ^ N , i n which the j-th particle has a constant drift fij, 



1 < j < N 19]. It is an extension of the Dyson model (JSJ) in the sense that, if we consider 
the scaled process eZ°(t/e 2 ),t > with e > for the case jl = and take the limit e — > 0, 
then the limit process is equivalent to X(£),£ > 0. He showed that the process Z^(£),£ > 



20 



22| 



is associated with the quantum Toda lattice with the Hamiltonian 

1 N f) 2 N ^ 

3=1 J 3=1 

The O'Connell process is very rich in mathematics connecting with quantum integrable 
systems, representation theory of Lie groups/algebras, the Whittaker functions, theory of 
intertwining relations of Markov processes, and so on. He discussed the importance of his 
process to study a model of 1+1 dimensional directed polymers in random environment with 



finite temperature 19 ]. 



The purpose of the present paper is to discuss the O'Connell process as a generalized 
version of vicious BM with appropriate conditions at least for the special case jl = 0. (See 



23 



for other generalizations and recent topics of vicious BM and noncolliding BM. 
We note that interesting connections between random growth models and the Toda lattice 
Hamiltonian is discussed in 



The paper is organized as follows. In Sec. II through the Feynman-Kac formula, we 
introduce a system of Brownian particles with the killing term which is in the same form 
as the potential term in the quantum Toda lattice Hamiltonian ([9]) and discuss it as a 
generalization of the vicious BM. In Sec. Ill the transition probability density Q7v(£,y|x) of 
the iV-particle system of killing BMs is expressed as an integral of a product of eigenfunctions 
of the quantum Toda lattice over the Sklyanin measure. Then asymptotics of Qnifi y| x ) m 
t — > oo is estimated (Lemma 1). In Sec. IV we introduce a drift ft in our iV-particle system 
of killing BMs and define the transition probability density of the killing BMs conditioned 
to survive up to time < T < oo. By taking the double limits T — > oo and jl — > 0, 
we obtain the transition probability density P/v(^,y| x ) for the killing BMs with /2 = 
conditioned to survive forever. The main theorem is given there (Theorem 2), by which the 



equivalence between the present conditional process and the O'Connell process with \x = 
is concluded. We discuss a one-dimensional diffusion process studied by Matsumoto and Yor 
in Sec.V as a motion of relative coordinate in the N = 2 case of our process. The 



Matsumoto- Yor process with \x = is realized as a one-dimensional killing BM conditioned 
to survive forever. In Sec. VI we discuss some distributions obtained by setting the special 
initial conditions. Section VII is devoted to summary and concluding remarks. Appendix 
lAlis given for proving an asymptotics used in Sec.V. Some details of the N = 2 case of the 
O'Connell process are given in Appendix IB1 



II. QUANTUM TODA LATTICE AND FEYNMAN-KAC FORMULA 

Let N G {2, 3, . . . }. Consider the eigenvalue problem of the quantum Toda lattice Hamil- 
tonian (jUJ), 

ftjv# 7 (x) = 7* 7 (x), x G R N . (10) 
For A = (Ai, . . . , Ajy) G C N , the eigenf unctions of (ITU1) with eigenvalues 



N 

1 - 

7 



*X>5 (id 



2 
j'=i 



have been extensively studied |20H22j . which are expressed by t/A^x) in the present paper. 
Let T denote a triangular array with size N, T = (Tkj, 1 < j < k < N). We consider that 
the N(N — l)/2 elements Tkj of T are independent variables and introduce a function of 
them as 

N f h fc-1 x 

fc=l Vi = l 3=1 y 

AT-1 fc 

_ | e _ ( r fc,i _T fe+ij) _|_ e -( T fc+i,j+i- T fe,j)| ; 



fe=l 3=1 



which depends on A = (Ai, . . . , A at). For a given x G K , let Ta^x) be the space of all real 
triangular arrays T with size N conditioned 

T N , j = x j , l<j<N. (13) 

We write the integral of a function / of T over T^x) as 

N k „oo N 

/ /(T)dT = Jin/ dT kJ f(T)l[5(T N/ -x e ). (14) 



Then the integral representation of i/j^ (x) is given by 

A 



^(x) = / (L) dT. (15) 

'rjv(x) 



This multivariate function is a version of Whittaker function (see [19[ and references therein). 

As a stochastic version of the Schrodinger equation of the quantum Toda lattice (obtained 
by performing the Wick rotation in the Schrodinger equation), we consider the following 
diffusion equation 

d 

—u(t,x) = £ N u(t,x) (16) 
with the infinitesimal generator of the process 

C N = -n N = ^A-V N (^), (17) 

where 

JV-l 

v N {*) = j2 e ~ (Xj+i ~ Xj) - ( l8 ) 

If we follow the method of separation of variables by setting u(t, x) = T(£)\l/ 7 (x), ( fTBl) is 
decomposed into the equations 

and (ITUj) . Then we can conclude that for any A G C N , 

N \ 

Vf>(x) (19) 

V" j=i / 
solves the diffusion equation (fl6l) . 

In the context of quantum mechanics, the function Vjv(x) given by (fl8l) plays, as a matter 
of course, a role of potential energy Then the quantum system prefers the state > Xj to 
the state < Xj, 1 < j < N — 1, since the former has lower energy than the latter. On the 
other hand, in the context of stochastic calculus, — Vjy(x) term in the infinitesimal generator 
of the process (TlTl) acts as a killing term. We consider N independent one-dimensional 
standard BMs starting from 0, Bj(t), 1 < j < N, and for x e R N set B x (t) = x + B(t), 
where each element B x - 3 (t) = + Bj{t) is a one-dimensional standard BM starting from 
2^,1 < j < ^V"- Then the Feynman-Kac formula (see, for instance, [ljj]) implies that the 
function 

Q^(t,y|x)=E l(B y (t) =x)exp|-^ lV(B y ( S ))^|J (20) 




solves the diffusion equation (fTBj) under the initial condition 



Qjv(0,y|x 



:) = *(x-y), 



(21) 



where E[-] denotes the expectation over all realizations of iV-dimensional Brownian paths, 
{B y (s) : < s < £}, starting from y, and is the indicator function of the event u; 

1(oj) — 1 if lo is satisfied, 1(lo) = otherwise. The function Q N {t,y\x) is the transition 
probability density of the process (Ti~6j) from a configuration x to a configuration y in time 
interval t > 0. In the Feynman-Kac formula (T2"U1) . we consider a collection of all paths of 
BM in Mr starting from y to x. (Though the time direction is backward, it is irrelevant 
in calculation, since BM is time-reversible.) The point of this formula is the following. In 
order to give the transition probability density Qjv(t,y|x), we have to put a weight 



to each realization of path of the N- dimensional BM and take a summation over all re- 
alizations of paths. It is obvious that takes a real value in [0,1]. Then this sum- 
mation of weighted paths (a path integral) can be identified with a statistical-ensemble 
average of Brownian paths, in which each path is included in the ensemble with proba- 
bility wn and is deleted with probability 1 — w^. Deletion of an iV-dimensional Brown- 
ian path is interpreted as an event that the N- dimensional BM is killed in the time in- 
terval [0, t]. The weight wn is then regarded as the probability that the particle in M N 
survives up to time t. (See Corollary 4.5 and explanation given below it in Chapter 4 
of [ijj for the equivalence of the Feynman-Kac formula to Brownian motion with killing 
of particles.) Eq. (l22|) gives the dependence of the survival probability on the realization 
of path {B y (s),0 < s < t}. If the iV-tuples of Brownian paths are "well-ordered" in 
the spatio-temporal plane, B\ 1 {s) < -Bf 2 (s) < •■■ < B^(s),0 < s < t, and moreover 
B V j+i{s) 3> B^ 3 (s), I < j < N — 1,0 < s < t, wn is large, while for a particle on the path 
{B y (s),0 < s < t} in which B^{s) < B y /(s), 1 < j < N for some s G [0,t], w N is small. 




(22) 



7 



If we introduce a parameter e > 0, then we can see that 

= l( Bf (s), . . . , B y N N (s) do not collide during [0, t] ) 

= 1 ^B y (s) G Wat, < Vs < t j . (23) 

In this sense, the process (fT6|) with (ITT)) and (|T8|) is an iV-particle system of killing BMs, 
which can be regarded as an extension of the absorbing BM in Wat. In the next section, we 
explain how to express the transition probability density given by the Feynman-Kac formula 
( l2"0j) as a superposition of the Toda lattice eigenfunctions (IT5|) . 

Remark. If we consider the present process not as an N- dimensional BM in M. N but as 
an iV-particle system of one-dimensional BMs, fl20|) gives the transition probability density 
in the case that mutual killing of particles does not occur at all in time duration t, since x 
and y are both iV-particle configurations, x, y G M. N . In order to discuss processes, in which 
mutual killing of particles actually occurs and total number of particles decreases in time, 
we have to specify the way how to choose pair of particles which are annihilated; e.g. the 
pair (J,j + 1) attaining min{5^Y(^) — B v k k {t)} is chosen. Note that in the original vicious 
BM, colliding pairs of particles are pair annihilated. In the present paper, however, we are 
interested in the process conditioned that all N particles survive. 



III. TRANSITION PROBABILITY DENSITY AND ITS LONG-TERM ASYMP- 
TOTICS 

The problem which is discussed here is how to determine the function g^(y) of A G 
C^, y G R N and a subset S of such that the integral of (TL9]) 

jf exp g Ajj (x)#(y)dX (24) 

is equal to Qn{^i y|x) given by fl20|) . This problem is solved by applying the theory of the 



Sklyanin measure 42| defined by 

(2vn 



s " {x)dx ^j^hw n {(A t -A j) sin ^' At) }n^ (25) 

' \<j<k<N ^ ' 1=1 



for A £ E 



where i 



-1. As a multi-dimensional extension of the fact that 



Macdonald's functions of imaginary order A £ zlR, Aj^e -1 ), make complete basis of a suitable 
set of functions with respect to the measure si(X)d\ = (i/ir) 2 \ sin(7rA)<iA (see p. 131 of {43) 
and Sec.V in the present paper); 



7T* 



/ioo 
K x (e- x )K x (e- y )\ sin(7rA)dA = 6{x - y), 
-ioo 



(26) 



x, y £ R, the following is valid [20|, |21| . 



4 7V) (x)^ ) (y)^(A)rfA = 5(x-y), 



(27) 



x, y £ K^. Note that if A £ (iR) N then ijj (1 V(y) is the complex conjugate of ip^^y)- 

— A A 

Since the transition probability density ( 1201 is a unique solution of the diffusion equation 
( TT5"j) satisfying the initial condition (12~T|) . the following is concluded; 



Qjv(*,y|x) 



*^^?/^w (x)v ,w (y)s7v(X)dX _ 



(28) 



It should be noted that Lemma 4.6 in [20j implies 



V^Wf } (x)rfx = * -1 £ %(A) - r) 



(29) 



for A, t £ (iM)^, where iSat is a set of all permutations of N indices and cr(A) = 
(Ao-(i), . . . , A ct (jv)) for a £ Sn- The orthogonality (129]) guarantees the Chapman-Kolmogorov 
equation for the transition probability density 



Q N (t, z\y)Q N (s, y|x)dy = Q N (t + s, z|x), x, z £ 



3 /v 



(30) 



for < s,t < oo. As a matter of course, (|20|) should satisfy it by the Markov property of 
BMs. 

If we change the integral variables in ( )28l) with ( 1251) as Aj i— >• za,- by Xj^t/2 = iuj, 1 < j < 
N, we obtain the following expression for Q^(t,y\x), 

2^/2 



Qjv(t,y|x 



(2^)^! 

w n 

l<;<fc<iV 



sinh{7rA/2/t(z/fc — Uj)} 



7TV2/t 



l\/2/tv —i\/2/tv 



■X N ) 



du, (31) 



9 



x, y G M , t > 0. This integral expression should be compared with the Schur function 
expansion of (jlj [lo| . 

t -N 2 /2 

q N (t,y\x) = ^ N/2 h N {x.)h N (y) 

x E fl (v . + N- jV exp (-^T^ ) ^(*MK(yM), (32) 

where {z/ = (z/ 1; u 2 , . . . ) : V\ > v 2 > . . . } are partitions of integers, £(z/) denotes the length 
of z/, and s u (x) is the Schur function 

We can prove the following asymptotics of the transition probability density 

Lemma 1 Let 

N 2 , , 

«jv = — (33) 



and C N = U" =1 T(n)/(2ir) N / 2 . Then 



for x,y G R N . 

Proof. In the limit t — > oo, (I3TJ) will behave as 

2iV 2 /2 



lim ^Q^t, y|x) = 4 N \x)4 N \y) (34) 



Qiv(^y|x)^-^r^f) (x) ^) (y) / e -M 2 JJ K -^.) 2 rfz7. (35) 

^ ^ ' ^ RJV l<j<k<N 

The integral is a version of the Selberg integral (a special case with 7 = 1 and a = 1 in Eq. 



(17.6.7) of [16|), 

TV 

/ e-l^l 2 J] (u k -u j f = (2n) N / 2 2' N2 / 2 l[nl (36) 

^ RJV 1<J<A,</V n=l 

Then the lemma is obtained. 1 

IV. O'CONNELL PROCESS AS AN iV-PARTICLE SYSTEM OF KILLING 
BROWNIAN MOTIONS CONDITIONED TO SURVIVE FOREVER 

Let \x = (fix, . . . , fi N ) G M. N and introduce a drift term in the diffusion equation 

J^(f,x) = /*ti*(f,x) (37) 



10 



with 



N 



dxj 

It is easy to confirm that if u(t,x) solves ffTET) . then 

vP{t, x) = exp f — -|/I| 2 + /2 • x j x) 



exp 



^ iV TV \ 

■ 2 E + E »i x i x ) ( 39 ) 
j'=i i=i / 



solves (j57|) . The formula (1551) is called the drift transformation from u(t, x) to u^(t,x) jl^ . 
Then we can prove that the transition probability density for the process with drift /I is 
given by 

Q&(t,y|x) = exp|-^|/2| 2 + /i-(x-y)|g J v(t,y|x) (40) 

for x,yelV>0,/Ie R*. 

ForxGl w ,0<T< oo, let 

A^(T,x)= / gJ(T, y|x)dy, (41) 
Jr n 

which gives the probability that all iV particles survive up to time T in the process of killing 
BMs with drift /X 

Now we consider the killing BM with drift fl conditioned to survive up to time T. Assume 
that < s < t < T and let Pj^ T (s, x; t, y) denote the transition probability density from 
a configuration x at time s to a configuration y at time t of this conditional process as 
illustrated by FigJH By the Markov property of BMs, we treat the paths only after time s 
for Pjy T (s, x; t, y). Since configurations at times s and t are both specified to be x and y, 
but configuration at the final time T is arbitrary, P£ T (s, x; t, y) should be proportional to a 
product of Q%{t — s, y|x) and Mj^(T — t, y). The product should be divided by J\fft(T — s, x) 
to give Pjy T (s, x; t, y), since J RN P^ T (s, x; t, y)dy = 1 should hold for any x G M. N . Here we 
can see that by the Chapman-Kolmogorov equation (I30l . f RN J\f£(T—t, y)Q^(t—s, y |x)dy = 
M^{T — s, x). That is, we have the formula 

Nft(T - s,x) 

11 



T 




A "J <* 
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FIG. 1: An illustration of the paths in the case N = 5 and /I = such that all particles survive 
up to time T, in which particle configurations at times s and t, < s < t < T, are specified to 
be x = (xi, X2, • • • , £5) and y = (2/1,1/2, • • • , 2/5), respectively. By the Markov property of BMs, the 
behavior of paths after time s is independent of that before time s. The total survival probability 
of the process during time interval [s, T] is Mj^(T — s, x) under the condition that the configuration 
at time s is x. If we fix the configuration at time t to be y, s < t < T, the survival probability 
density is given by a product M^{T — t, y)Q%(t — s, y|x), since the events before time t and after 
time t are independent by the Markov property of BMs. As shown by (142 j) . the ratio of them gives 
the transition probability density P& T (s,x;t,y). 

By definition, this conditional process is a temporally inhomogeneous process, which is 
clarified by the fact that the transition probability density (1421) is the function not only of 
the time difference t — s but also of T — s and T — t. 

In order to obtain the temporally homogeneous process, we take the limit T — > 00. 
Here we note that ipo N \y) is the eigenfunction of the quantum Toda lattice (Q with zero 
eigenvalue. If y J+ i ^> yj for all 1 < j < N — 1, the potential energy f|T8|) becomes zero and 
then (y) behaves similar to the harmonic function /ijv(y). On the other hand, it is known 



12 



ihat ipf\y) ~ exp{-2e-fe+ 1 -w ) } ->• as yj+i-yj -»■ -oo, 1 < j < iV-1 for any A G 
2o| . Then if /I is chosen so that jl G Wjv, then the integral / Rjv ipQ^iy) exp(— /2 • y)dy is 
finite. In such a case, by Lemma 1 and the drift transformation fj40l) . 

r r a N 

lim el^A/^T - s, x) = exp(/2 • x^x^iV), (43) 

T->oo Otv 



< Vs < oo, with 



lg(N)= f 4 iV) (y)exp(-//.y)dy, /2 G W^v- (44) 



Then we have 



P#(*-s,y|x) = lim P^ T (s,x;t,y) 

T— i>oo 



exp{-/2 • (x - y)} ° {N) , ^ Qn( 1 - s , yl x 



= eH^&Q w(i - S ,y|x) (45) 

YO ( X ) 

for < s < t < oo, x, y G , jl G Wjy. 

Finally we take the limit jl — > 0,/2 G Wat- The obtained transition probability density is 
given by 

P N (t,y\x) = lim P#(t,y|x) 

M^o,/JeW]v 

= WT^(^I X )' ( 46 ) 
( x ) 

x, y G M^, t > 0. It is an extension of fl6]), where the Vandermonde determinant fojv(x) 
is replaced by the eigenfunction i/jq N \x) with zero eigenvalue of the quantum Toda lattice. 
Now we state the main theorem of the present paper. 

Theorem 2 The function (|46|) . which is obtained as the transition probability density of 
the killing BM conditioned to survive forever, solves the following differential equation, 

^u(t, x) = ~ Au(t, x) + V log 4 N) (x) • Vu(t, x) 

3=1 J J=l J J 

under the initial condition 

u(0,x) = 5(x-y). (48) 
13 



Proof. By (06 



9 p (t vixi - ^° (y) ^(x) g , 4 N \y) d 



dx) (4 } (x)) 3 V dx * ) 

(^(x)) 2 9x] WM ' y| j ( ^) (x))2 

, } (y) 5 2 



Since 



/ (at) / .njQ^yW- ( 49 ) 

9 log 4^ (x) i #r } ( x ) 



9xj 4 N \x) dx i ' 

RHS of (HU) is equal to 

m ( x )r 3=1 °^ V^o (x) 2 j= i ax i 

Since ^ (x) is the eigenfunction of jHJ) with zero eigenvalue, 

i # t& iV*) ( \ at-1 

^W(x) = -i ^ ^« + £ >-*"^P(»l = o. (si) 

i=i J 3=1 

and Qjv(t,y|x) satisfies the equation 

-g w (t, y |x) = - ]T ^ Q*(*, y|x) - X) c-^-^O^, y|x), (52) 

3=1 Xj 3=1 

we can see that (1501 is equal to dPjy(t, y|x)/<9t; that is, Pjv(^y| x ) satisfies ( HT)) . Since 
Qjv(0, y|x) = <5(x — y), (T4"6]) gives P/v(0, y|x) = S(x — y). The proof is then completed. | 



19] introduced a diffusion process in M. with the infinitesimal generator of 



O'Connell 
the process 

^S = iA + Vtog4 JV) (x).V (53) 

with /2 G R . Our theorem states that the special case with /2 = can be realized as 
a vicious BM, which is obtained above as a system of killing BMs conditioned to survive 
forever. Corresponding to the backward Kolmogorov equation (j47p . the SDEs of the process 
Z(t) = (Zx{t), . . . , Z N (t)),t > is then given by 

dZ^i) = dBj(t) + Fj N) (Z(t))dt, 1 < j < N, t > (54) 

14 



with 



d log tp t 



(AO/ 



(55) 



where Bj(t), 1 < j < AT, t > are independent one- dimensional BMs. 



V. MATSUMOTO-YOR PROCESS AS A KILLING BROWNIAN MOTION CON- 
DITIONED TO SURVIVE FOREVER 



For N = 2, the eigenfunction (1T5|) of the quantum Toda lattice is 



JT 2 (x 1 ,x 2 ) 



r(2) 



exp 



AiTx,! + A 2 (x! + x 2 - T ltl ) - {e-C^- 3 *) + e^ 2 -^ 1 )}]^ 

dTx,!, (56) 



e A 2 (, 1+ , 2) / e (Ai-A2)ri,iexp 



+ 



(Ai,A 2 ) G C 2 ,(xi,x 2 ) G M 2 . Change the integral variables T^i n- s by e TM = e^ 1+X2 ^ 2 s. 
Then 



WaU)(*i^2) = e^+^X^V 2 / s^-^^exp 



. e -(* a -xi)/2 



tZs. (57) 



Let I,,(z) be the modified Bessel function of the first kind 



z/2) 



fc=0 



(fc + i)r(fc + z/ + i) 



, |z| < oo, |argxr| < 7r 



(58) 



and K„(z) be Macdonald's function 



43] 



Kv {z) = l 1 -^}, (a^:| < -. fur v + 0. ± 1. ±2. 



2 sin(z/7r) 



(59) 



and for integers v — n, 



K n (z) = \im K v (z), n = 0,±l,±2, 



(60) 



which are both analytic functions of z for all z in the complex plane cut along the negative 
real axis, and entire functions of v. We can see that I u {z) and K u (z) are linearly independent 
solutions of the differential equation 



d 2 w 1 dw 



dz 2 z dz 



v 



+ -— - 1 + — tw = 



(61) 
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For x > and v > 0, /^(x) is a positive function which increases monotonically as x — >■ oo, 
while K v (x) is a positive function which decreases monotonically as x — > oo. Since K v (z) 
has the following integral representation 



KJz) 



s u 1 exp 



z / 1 

2 S + J 



(is, 



( 1571) is written as 



(62) 



(63) 



The infinitesimal generator of the process fl53l) is then given for iV = 2 as 



»(M1i/42) 



| 2 ^ 

9 A + E^ 10 ^Sa 2 (^ ; 



«9 2 <9 2 
+ 



^2, 



2 V&r 2 



(9a; 2 



1, , / 9 9 

+ 2 (/il + /i2) fe + 



-(x 2 -a;i)/2N 



-(x 2 -xi)/2 



(9a;2 
8 



Kn-npe-l**-*)/ 2 ) \d Xl dx 2 , 

where K' v {z) = dK u (z)/dz. If we change the variables (xi,x 2 ) (£,rj) by £ 



(64) 



v = - 

where 



(xx + x 2 )/2, 

[X! - x 2 )/2 - log 2, (ED) is decomposed into two parts, g^>i*) = (^w+w + Q^^)/2, 



r u _ Id 2 d_ 

y ° ~ 2de dC 



1 d 2 



MY 



Id 2 d n TJr . d 



d 
drj 



(65) 
(66) 



The former is the infinitesimal generator of the one-dimensional BM with a constant drift 
v = \i\ + fi 2 and the latter is that of the diffusion process studied by Matsumoto and Yor 



with parameter \i — \i\ — fi 2 40|, |41| . It implies that, in the N = 2 case of the O'Connell 



process with parameter /2 = (fii, fi 2 ) e R 2 , the center of mass £ behaves as (a time change 
1 1 — y 2t of) a BM with a drift \i\ + \i 2 and the relative coordinate rj behaves as (a time change 
1 1 — y 2t of) the Matsumoto- Yor process with parameter n\ — fi 2 . 



For fj,eR, let 



with 



1<P_ 
2dx^ 

V(x) 



V(x) — \i 



d 

dx 



-2x 



(67) 



(68) 
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It is the infinitesimal generator of the one- dimensional drifted BM with a killing term —V{x). 
The transition probability density is given by 

Q*(t,y\x) = E l(B™(t)=x)exp{~J e~ 2B "' y(s) dsj 

= e -^t/2+^-y) / e ^K x (e- x )K x (e- y )Xsm(7cX)dX, (69) 

7T / 

" J —ioo 

where B^' y (t) = y + B(t) + fit with the one-dimensional standard BM, B(t), starting from 
0; B(0) = 0. 

The survival probability up to time T > of this one-dimensional killing BM with drift 
fi is given by 

/oo 
Q»(T,y\x)dy (70) 
-oo 

for the initial position We can prove that, if fi > 0, it has the long-term asymptotics, 

lim A /^T 3/2 e^ T/2 ^(T,a;) = e^K^) [°° ^(e^e^dy 

V 2 J-oo 

= 2^- 2 (T(fi/2)) 2 e^Ko(e- y ). (71) 

The proof of (ITTj) is given in Appendix |A] Then the transition probability density for this 
one- dimensional killing BM conditioned to survive forever is given in the limit fi — > 0, fi > 
as 

A/WT — t v) 
P(t,y\x) = lim lim J \) . \ y) Q»(t,y\x) 



K (e 



/ioo 
e x2t/2 K x (e- x )K x (e- y )\ sin(vrA)rfA, (72) 
•ioo 



K (e x ) n 

x, y G K, t > 0. It is easy to see that (I72|) satisfies the diffusion equation 

= G^ Y u(t,x), t>0, (73) 

under the initial condition u(0, x) = 5{x — y). 

Matsumoto and Yor showed that the stochastic process 

Z^ Y {t)= log y\ 2B ^ds\ - B»(t), t>0 (74) 

with B^(t) = B^'°(t) = B(t)+fit is a diffusion process, whose infinitesimal generator is given 
by ( 1661) for any /iGE. Here we have shown that, when fi = 0, the Matsumoto- Yor process 
( 1741 can be constructed as a one-dimensional killing BM conditioned to survive forever. 

See Appendix B for more detail on the relation between the Matsumoto- Yor process and 
the N = 2 case of the O'Connell process. 
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VI. SPECIAL INITIAL CONDITIONS 



In this section, first we consider the one-dimensional diffusion process with the infinites- 
imal generator ([6"7]) with ( IBS]) . Let < T < oo. Then the transition probability density of 
the process conditioned to survive up to time T is given by 

~ ( IZoQ fl ( T ~~ s,z\x)dz' (75) 

< s < t < T, x, y G M., where and A^* are given by ( 169]) and ( 1701 . respectively. The 
asymptotics of K\(e~ x ) in x — > — oo is independent of A |43| 



K x (e x ) ~ y ^^exp(-e x ) as x -> -oo. (76) 

Then if we define 

P^,y|-oo)= lim P£(0,x;t,y), (77) 

it is given by 

P£(t,y\ -oo) = 2 ,- 2(r ^ /2))2 T3/2 ^ (T " t)/2 ^ 2/T ^-4t)e-^A-(T - f, y ), (78) 



where r (t) is given by (|A9]) and (1A10I) in Appendix A |44|. Since A/" M (0, y) = 1 by definition, 
we obtain the distribution at time t = T, 

P£(T,y\-oo) = c >1 (T)6 e ~y(T)e- fiy J y G E, T>0 (79) 

with c^(T) = V2^2 2 -^(T(fi/2))- 2 T 3 / 2 e- n2 / T . On the other hand, by (171]) . if we take the 
temporally homogeneous limit T — )■ oo in (178]) . then we obtain the distribution 



P l *{t,y\ -oo) = lim P£(t,j/| -oo 



T^oo 



= 2e- /A/2 fl e - ! ,(t)fro(e- 3 0, !/^> * > 0. (80) 

Next we consider the A-particle system of the killing BMs with drift jl G Wjy conditioned 
to survive up to time T, < T < oo. The transition probability density is given by ( 142]) . 
Corresponding to (1761) . the asymptotics of (x) in Xj — >■ — oo, 1 < Vj < A is independent 



of A (see Remark 8.1 in 



19j). Then we obtain 



p 5,t(*> y| - °°) = PjL^x^y) 

oo,l<j<JV 



e lM| 2 (T-t)/2 

J*(A T) 9jv(t ' y) GXP(-/2 ' y)A/ ^ (T ~ £ ' y) ' (81) 
18 



where 



with 



and 



ji*(N,T)= f e T ^ x V 2 s N (X)I^JN)d\ (82) 
I> 1 JN)=[ ^ ( ^ ) (z)exp(-/i-z)ciz, Xe(iR) N , fIeW N , (83) 

~ A Jrn " A 



e^(t,y)= / e t ^ x V 2 ^(y)s N (X)dX. (84) 

J(iM) N 

At time t = T, ( 18~T|) gives the distribution 

P% )T (T, y | - oo) = * e N (T, y) exp(-/i • y) , y6l ff , T > 0. (85) 



On the other hand, if we take the limit T — >• oo, fjSlH gives the distribution 



Py(t,y| - oo) = Jim P^ T (t,y| - oo) 

= e-i^/ 2 e JV (t, y )4 7V) (y) ) yeK ff , t>0. (86) 



The three-dimensional Bessel process is defined as the radial part of the three-dimensional 
BM and abbreviated to BES(3). The above results will be compared with the Imhof rela- 



S(3) and the process called a meander and its multivariate generalizations 



tion between B 
discussed in [ol 1(3] . 



VII. SUMMARY AND CONCLUDING REMARKS 

The vicious BM is obtained as a diffusion scaling limit of Fisher's vicious walk model 

nag 

. It is an iV-particle system of BMs in one dimension, whose positions are arranged 
in the order x± < x 2 < ■ ■ ■ < x^, such that if and only if two neighboring Brownian particles 
collide with each other then they are pair annihilated, while they can enjoy free Brownian 
motions if they are all located separately from each other. In the present paper we have 
considered a system of iV Brownian particles with the killing term 

JV-l 

- y JV ( x ) = - e-to+ 1 - x i\ (87) 
3=1 

That is, the interactions between neighboring Brownian particles are long-ranged and the 
risk to be pair annihilated exists always, which is expressed by a rapid decreasing function 
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(1ST)) of the distance of the two particles Xj+i —Xj. We regard this system of mutually killing 
BMs as a generalized version of vicious BM, since the original vicious BM can be identified 
with the system of BMs with the killing term obtained by — lim e ^ ,e>o ^/v( x /£)- 

Though the original vicious BM has only contact interactions, if we consider the system 
conditioned never to collide with each other, then we obtain a system of BMs with long- 
ranged interactions; the SDE is given by Eq.flSJ), in which between any pair of particles 
there acts a repulsive force proportional to the inverse of distance of the pair js|. This N- 
particle process is equivalent to the eigenvalue process of an N x N Hermitian-matrix- valued 
BM introduced by Dyson in order to dynamically simulate the eigenvalue statistics of the 
Gaussian unitary ensemble (GUE) of random matrices (the Dyson model) 15 l^r| . 



As discussed in 



45 



46], the equivalence between the eigenvalue process of Dyson and the 
noncolliding BM (the vicious BM conditioned never to collide) is the iV-variate extension of 
the equivalence between BES(3) and the one-dimensional BM conditioned to stay positive. 
In this sense, the Dyson model can be regarded as a many-particle generalization of BES(3). 

Apart from the equivalence between the BES(3) and the conditional BM to stay positive, 
the following equivalence is established. Let M(t) = maxo< s <t B(s), t > 0, and define a 



process Y(t) = 2M(t) - B(t),t > 0. Then Y 



as Pitman's '2M — X' theorem 



t) is equivalent to BES(3), which is known 



Pitman's '2M — X' theorem, anot 



BM) has been reported 25 



26 



47) | (see also 40j, |4l|, |48j). As a multivariate extension of 
rer construction of the Dyson model (the noncolliding 



28 



31 



32] 



Matsumoto and Yor studied the stochastic process Z!^ Y (t),t > given by ( ]T4"j) . We can 
see that 



lim eZ^ v {t/e 2 ) 



£— >-0,£>0 

= lim 

e^0,£>0 

= 2max5 /1 (s) - B»(t), t > 0. 

0<s<* 



e\og < 









(88) 



Then, when fi = 0, this e — > limit is equivalent to Y(t) and thus with the BES(3). In this 
sense, the Matsumoto- Yor process is a generalization of the BES(3) 40, |41| . 

O'Connell Q introduced an X-particle process Z^(t) = (Z?(t), Z%(t)),t > 0,/2 G R N 
as a multi-dimensional generalization of the Matsumoto- Yor process. Corresponding to the 
fact that the Matsumoto- Yor process is a generalization of the BES(3), the O'Connell process 
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is a generalization of the Dyson model. Actually, he showed that \im £ ^ 0t£> oeZ°(t/e 2 ), t > 
is equivalent to the Dyson model in the sense of an extension of Pitman's '2M — V theorem 



191 ] . We pointed out that the BES(3), the original vicious BM, and the Dyson model (the 



noncolliding BM) can be regarded as ultradiscretizations 49j of the Matsumoto-Yor process, 
the BMs with the killing term in the same form as the quantum Toda lattice potential, and 
the O'Connell process, respectively 

In the present paper, we discussed another construction of the O'Connell process apart 
from the extension of Pitman's '2M — V theorem. In the special case with fl = 0, we have 
shown here that his process is given as a generalized version of vicious BM conditioned to 
survive forever. In order to demonstrate that the relation between the present generalized 
vicious BM and the O'Connell process is a multivariate generalization of the relation between 
a killing BM and the Matsumoto-Yor process, we showed in Sec.V that the Matsumoto-Yor 
process with /i = is obtained as a killing BM conditioned to survive forever. 

We want to emphasize that the present analysis is indeed based on the idea of O'Connell 
to discuss interacting diffusive particle systems using the exact solutions of the quantum 
Toda lattices |l^ |. 

In Sec. I in the present paper, we listed up three fundamental properties of the noncol- 
liding BM. They are all inherited by the O'Connell process in the extended form, (i) The 
Karlin-McGregor determinantal expression of qx{t, y| x ) ; which is expanded by the Schur 
functions (l32l) . is generalized by the integral formula (l3Tl) . (ii) The harmonic transform I3I ] 
from to pn (EJ) by the harmonic function hjsr given by the product of differences of vari- 
ables (the Vandermonde determinant) (j2J) is now given by the formula (1461) from Qn to 
P/v- There h N is replaced by an eigenf unction of the infinitesimal generator £ N (the 



Hamiltonian Hn of the quantum Toda lattice). (See |50M52l| for harmonic transforms of one 
dimensional generalized diffusion processes.) (iii) Theorem 2 in Sec. IV gives the extended 
version of the Kolmogorov equation of (El). 

There are a lot of future problems. In noncolliding diffusion processes, if we study the 
situations starting from "the all zero state" and observe particle distributions at an arbitrary 
time < t < 00 in temporally homogeneous processes, and at the ending time t = T in 
temporally inhomogeneous processes defined only in an finite time- interval [0, T], we have 
obtained the eigenvalue distributions of random matrices in a variety of ensembles studied 
in random matrix theory 10] . In Sec. VI, we demonstrated that "the all —00 state" and 
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temporally inhomogeneous versions of processes will play important roles in the O'Connell 
process. The noncolliding diffusion process is determinantal, in the sense that for any finite 
initial configuration all multitime correlation functions are given by determinants associate 
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541. It will be a challenging problem 
1171 ] and solvability are 



with an integral kernel called the correlation kernel 
to clarify how matrix-structures (i.e. symmetries of systems 
inherited by the family of O'Connell processes. 
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Appendix A: Proof of (I71|) 

Let Jq(z) be the Bessel function of the first kind of order 0, 

(-l) k (z/2) 2k 

Since the equality 



J o(*) = E v > N<°°- ( A1 ) 



7T f°° 

K x (x)K x (y) = — — — / J {y/2xy coth u - x 2 - y 2 ) smh(uX)du (A2) 
2 sm(vrA) Ji os ( y / x ) 



og{y/x) 

holds for x > 0,y > 0, |ReA| < 1/4 [43], ([69]) is written as 



poo 

Q^{t,y\x) = e -M 2 */2+M(^) / du j (y/2e-( x +v) coshw - e~ 2x - e~ 2 y) 

^ J x-y 

/ioo 
dX e A2 * /2 Asinh(uA) 
-ioo 

POO 

/ uJ (y/2e-( x +v) cosh u - e- 2x - e' 2 " )e~ u2/2t du, (A3) 

J x—y 



' —ioo 
-/i 2 t/2+fi(x—y) poo 



27Tt 3 / 2 
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where we have performed the integral of A over zR. This gives 

,,y\x) 



lim ^t^e^Q^t, 
;— >oo y 2 

/ uJoC^e-^+f) coshw - e" 2 * - e" 2 ^. (A4) 
2 J x—it 



e 



-v 

We find that the A — > limit of ( 1A2I) gives the equality 

1 /-do 

iT (^)^o(y) = « / Mj (A/2e-( :r +2') cosh u — e~ 2x - e- 2y )du, (A5) 

and then (1A4I) gives 

lim M^e^Q^t, y\x) = e^ x '^K ( e - x )K (e- y ). (A6) 
t— >oo y 2 



The asymptotics of -Ko(e y ) is known as 43| 



r . _„. , log(2/e y )~y asy^oo 

^o(e y ) ~ <( _ (A7) 

/7r/ (2e~ y ) exp(-e y ) -)■ as y -)■ -oo. 

Then if yu > 0, K (e~ y )e~^ y dy < oo. Actually, for // > 0, this integral is the Mellin 
transformation of Kq(z) and we obtain 

/oo poo 
K {e- y )e~™dy = / K (z)z^ 1 dz = 2^ 2 (r(/j/2)) 2 . (A8) 
oo JO 

Therefore ( J7TT) is valid. ■ 

We note that, for the function 

• pioo 

e r (t) = — e xt/2 K x (r)\sm(7r\)d\, r > 0, (A9) 
27T^ J_ ioo 



Yor gave the following expression (see Eq.(6.b") on page 43 of [44j) 



r {t) = {2 J H)l/2 e n2/2t e~^ /2t e- rcosh ^smh rj) sin (^) dry, r > 0. (A10) 
Using this expression, Matsumoto and Yor reported the asymptotics (see Eq.(2.11) in 41]), 



t— >oo 



lim V27r^0 r (t) = K (r), r > 0. (All) 

Since the equality 



oc 



Q"(t,y|aO = e~^ 2+ ^ I exp <j -t - l(e" 2 * + e"*)} 8 e - (x+y)/s (t)- (A12) 



o 



2 2s v y J y s 

is established, the limit ( 1A6|) can be concluded also from flAllj) 
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Appendix B: N = 2 case of the O'Connell process 

By the equations (1251) . (128|) and (|63|) . we obtain 



2 , \2\ 



»oo »/ — lOO 



g 2 (t, y|x) = -Ar / dA x / rfA 2 e^ +A 2) t / 2 e ( Al+A2 )« :Cl+:!;2 )-^ + ^»/ 2 



X 



X Al _ A2 (2e-^-^)/ 2 )ir Al _ A2 (2e-fe^)/ 2 )(A 1 - A 2 ) sin^ - A 2 )}, (Bl) 



x, y G M, 2 ,t > 0. If we change the integral variables (Ai,A 2 ) (A, v) by A = Ai — A 2 , v = 
Ai + A 2 , we can calculate the integral with respect to v. The result is expressed by using 
the transition probability density of the Matsumoto-Yor process fl69|) with \i = as 

Q 2 (t, y |x) = p(2t, Vl + y 2 \ Xl + x 2 )Q° Q, - log 2| - log 2^ , (B2) 

where p(t,y\x) = e~^ y ~ x ^ 2 ^ 2t / \f2ixt. Therefore, from the long-term asymptotics ( 1A6I) of 
we can obtain the long-term asymptotics of Q 2 as 

Q 2 (i,y|x) ~ ^K (2e-^- x ^ 2 )K (2e-^- x ^ 2 ) 

7Tt 

= ^4 2, (x)4 2) (y) ast^oo, (B3) 
which coincides with the N = 2 case of (1341) in Lemma 1. 
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